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Abstract. The diffraction image of a quasicrystal admits a finite group G as a 
symmetry group, and the quasicrystal can be regarded as a quasiperiodic packing of 
copies of a G-cluster C, joined by glue atoms. The physical space E containing C can 
be embedded into a higher-dimensional space R fc such that, up to an inflation factor, C 
is the orthogonal projection of the set {(±1, 0, 0), (0, ±1, 0, 0), ... (0, 0, ±1)}. 
The projections of the points of Z k lying in the strip S = E+[— 1/2, l/2] k +t obtained 
by shifting a hypercube [—1/2, l/2] k +t along E is a quasiperiodic packing of partially 
occupied copies of C, but unfortunately, the occupation of clusters is very low. In our 
modified strip projection method we firstly determine for each point x € Z fe n S the 
number n(x) of all the arithmetic neighbours of x lying in the strip S, and project the 
points of Z fe n S on E in the decreasing order of the occupation number n(x). In the 
case when n(x) represents more than p% of all the points of the cluster C we project 
all the arithmetic neighbours of x (lying inside or outside S) . We choose p such that 
to avoid the superposition of the fully occupied clusters. The projection of a point x 
with n(x) less than p% of all the points of the cluster C is added to the pattern only 
if it is not too close to the already obtained points. 
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1. Introduction 

Quasicrystals are materials with perfect long-range order, but with no three-dimensional 
translational periodicity. The discovery of these solids in the early 1980's and the 
challenge to describe their structure led to a great interest in quasiperiodic sets of 
points fUEJ- The diffraction image of a quasicrystal contains a set of sharp Bragg 
peaks invariant under a finite non-crystallographic group of symmetries G, called the 
symmetry group of quasicrystal (in reciprocal space). In the case of quasicrystals with 
no translational periodicity this group is the icosahedral group Y and in the case of 
quasicrystals periodic along one direction (two-dimensional quasicrystals) G is one of 
the cyclic groups C 8 (octagonal quasicrystals), C w (decagonal quasicrystals) and C% 2 
(dodecagonal quasicrystals). 

Real structure information obtained by high resolution transmission electron 
microscopy suggests us that a quasicrystal with symmetry group G can be regarded 
as a quasiperiodic packing of copies of a well-defined G-invariant cluster C, joined by 
glue atoms |7j. From a mathematical point of view, a G-cluster is a finite union of orbits 
of G, in a fixed linear representation of G. A mathematical algorithm for generating 
quasiperiodic packings of interpenetrating copies of G-clusters was obtained by author 
in collaboration with Verger-Gaugry several years ago pQ . This algorithm based on strip 
projection method [01 El El HH] works for any finite group G and any G-cluster, but in 
the case of a multi-shell cluster the dimension of the involved superspace is rather high, 
and the occupation of the clusters occurring in the generated pattern is too low. Some 
mathematical results recently obtained by author [31 [S] simplify the computer program 
and allow to use strip projection method in the superspaces required by this approach. 
Now, our aim is to present a way to increase the occupation of clusters occurring in the 
generated quasiperiodic set. 

2. Two-dimensional packings of clusters 

Let G be one of the cyclic groups C$, Cw, C\ 2 . Each group C n can be defined as 

C n = ( a | a n = e ) = { e, a, a 2 , a n_1 } (1) 

and the formula 

/ f 2ir n . 2n . 2tt n 2n\ 

a{a, p) — [a cos p sin — , asm hp cos — (2) 

\ n n n n J 

define an M-irreducible representation in M 2 . The orbit generated by (a, (3) ^ (0,0) 

C n {a,/3) = {(a,f3),a(a,(3), a 2 (a,(3), a n -\a,f3)} (3) 

contains n points (vertices of a regular polygon with n sides). Let 

C 2 = {vx, v 2 , v k , -vi, -v 2 , -v k } (4) 
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where v x = (v n ,v 2 i), v 2 = (vi 2 ,v 22 ),..., v k = (vx k ,v 2k ), be a fixed G-cluster, that is, 
a finite union of orbits of G. From the general theory pQ (a direct verification is also 
possible) it follows that the vectors 

Wi = (^11,^12, ■■■,v lk ) and w 2 = {v 21 , v 22 , v 2k ) (5) 

from R k are orthogonal and have the same norm 

(wt, w 2 ) = VnV 2 i + v 12 v 22 + ■■■+ v lk v 2k = 

(6) 



iFill = V v n + v i2 + - + v ik = \\ w 2\\- 
We identify the physical space with the two-dimensional subspace 

E 2 = { awi + (3w 2 | a, f3 G R } (7) 
of the superspace R k and denote by E^ the orthogonal complement 

Ej = { i G R fc (x,y) = for all y G E 2 }. (8) 
The orthogonal projection on E 2 of a vector x G R k is the vector 

« 2X= U^)^ + U^)^ (9) 

\ K I K \ K I K 

where k = \ \wi\ \ = \\w 2 \\, and the orthogonal projector corresponding to E^ is 

tt^ : R k — > E^ vr^x = x-tt 2 x. (10) 

We describe E 2 by using the orthogonal basis {k~ 2 wi, k~ 2 w 2 }. Therefore, in view of 
(0) the expression in coordinates of 7r 2 is 

n 2 :R k — > R 2 n 2 x=((x,w 1 ),(x,w 2 )). (11) 

The projection W 2j fc = ^(A*.) of the unit hypercube 

A fc = |(xi,x 2 , ...,x k ) -^<Xi<^ for all % G {1,2, k} J (12) 

is a polyhedron (called the window of selection) in the (fc— 2)-dimensional subspace E^, 
and each (k — 3)-dimensional face of W 2j fc is the projection of a (/c — 3)-dimensional face 
of A fe . The vectors e x = (1, 0, 0, 0), e 2 = (0, 1, 0, 0), e k = (0, 0, 0, 1) from R fc 
form the canonical basis of R k , and each (k — 3)-face of A k is parallel to (k — 3) of these 
vectors and orthogonal to three of them. There exist eight (k— 3)-faces of A k orthogonal 
to the distinct vectors e^, ej 2 , ej 3 , and the set 



x = (xi,x 2 , ...,x k ) 



(13) 



Xi G {-1/2, 1/2} if i G {h, i 2 , i 3 } 
Xi = if % G" {ii, i 2 , i 3 } 

contains one and only one point from each of them. There are 

k(k-l)(k-2) 

6 1 ' 

sets of 2 3 parallel (k— 3)-faces of A k , and we label them by using the elements of the set 

1 2>k = {(ii,i 2 ,i 3 ) e Z 3 | 1 < h < k-2, h + l<i 2 <k-l, i 2 + 1 < «s < }. (15) 
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In R 3 the cross-product of two vectors v = (v x ,v y ,v z ) and w = (w x ,w y ,w z ) is a 
vector orthogonal to v and w, and can be obtained by expanding the formal determinant 

k 

(16) 



V X w 



ID, 



Wy W Z 



where {i, j, k} is the canonical basis of M 3 . For any vector u = (u x ,u y ,u z ), the scalar 
product of u and v x w is 



u(v x w) 



u. 



w„ 



u„ 



V 



y 



Wy UJ Z 

In a very similar way a vector y orthogonal to I — k — 1 vectors 
Ui = (un,Ui2, u i3 , ...,u ik ) i £ {1,2,3,...,/} 
from IR fc can be obtained by expanding the formal determinant 
d e 2 e 3 ... e k 

Un Uu Ui 3 ... Uik 
y = U 2 l U 2 2 U 2 3 ... U 2 k 

Ull U i2 Ui 3 ... Ui k 

containing the vectors of the canonical basis in the first row. 
(xi, X2, Xk) £ R fc , the scalar product of x and y is 



w, 



(17) 



(18) 



(19) 



For any x 



Xi 


x 2 


X 3 . 


■ x k 


Un 


Ul2 


Ul 3 ■ 


■ Uik 


«21 


U22 


U 2 3 ■ 


■ U 2 k 


un 


Ul2 


Ul3 ■ 


■ Uik 



(20) 
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(21) 



is a vector orthogonal to the vectors e^, e$, ek, Wi, w 2 , and 
(x,y) = (-l) k - 1 



Xi X 2 X3 
Vu V12 V13 
V21 V 2 2 V 2 3 



(22) 
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Figure 1. Left: The strip Si_3 and the window Wi^ in the case of a ID physical space 
Ei embedded into a three-dimensional superspace. Centre: A one-shell Cs-cluster Ci. 
Right: A fragment of a set f2 defined by C2. The nearest neighbours of any point q of 
fl belong to q + C%, which is a copy of C2 with the center at point q. 



for any x G M fc . The vector y belongs to E^, and since — it 2 ei is a linear combination 
of Wi and w 2 , it is also orthogonal to n^e^, 7r^e 5 , 7i 2 ek- Therefore, y is orthogonal 
to the (k — 3)-faces of W 2)fc labelled by (1, 2, 3). Similar results can be obtained for any 

(21,22,23) e J 2 , fc . 

Consider the strip corresponding to W 2jfc (see figure 1) 



>2,fc 



{a; G 



7T 



^ e w 2ifc } 



and define for each (21,22,23) G X 2j fc the number 



di x i 2 i z 



max 

aj £{-l/2, 1/2} 



«1 


"2 


« 3 
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(23) 



(24) 



A point x G M fc belongs to the strip S 2j fc if and only if 
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— ^n*2«3 
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(21,22,23) G X 2ifc 
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The set defined in terms of the strip projection method [HI El El E]] 
n = vr 2 (S 2ifc n Z fc ) = { vr 2 a; I 16 s 2 , fc n z fc } 



(26) 



can be regarded as a packing of translated partially occupied copies of C 2 . Since 

7T2ej = ((ei, wi), (ei, w 2 )) = (vu, v 2i ) = v { (27) 

we get 

7r 2 ({ x ± e x , x ± e 2 , sc ± e fc } n S 2jfc ) 

C {7T 2 X ± Ui, TT 2 X ± f 2 7T 2 X ± Wfc} = 7T 2 X + C 2 (28) 

that is, the neighbours of any point tt 2 xeQ belong to the translated copy ir 2 x+C 2 of C 2 . 

A larger class of aperiodic pattens can be obtained by translating the strip S 2j fc. 
For each t G M fc the set 

Q = 7T 2 ((t + S 2>k ) nZ k ) = { n 2 x I x-t G S 2>fc and x G Z k } (29) 
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Figure 2. Left: A one-shell Cio-cluster and a fragment of the corresponding 
quasiperiodic set. Right: A fragment of the quasipcriodic set defined by a two-shell 
Cio-cluster, obtained by using strip projection method in a ten-dimensional supcrspace. 
The starting cluster is a covering cluster, but for most of the points the occupation is 
extremely low. 

is a packings of partially occupied copies of C 2 . Some particular examples can be seen 
in figures 1-3. 

3. Modified strip projection method 

The algorithm based on the strip projection method presented in the previous section is 
very efficient. Hundreds of points of Q can be obtained in only a few minutes for rather 
complicated G-clusters C 2 , but, as one can remark in figures 1-3, for most of the points of 
Q the occupation of the corresponding cluster is very low. On the other hand, the images 
concerning the quasicrystal structure obtained by high resolution transmission electron 
microscopy show the presence of a significant percentage of fully occupied clusters. 

The number n(x) of the neighbours of a point n 2 x G f2 occurring in Q corresponds 
to the number of the points of the set {x ± e±, x ± e 2 , x ± e^} belonging to the strip 
t + S 2 ,fc. In our modified strip projection method we firstly determine n(x) for all the 
points of Z fe lying in the fragment of the strip we intend to project, and then we project 
the points in the decreasing order of the occupation number n(x). In the case when n(x) 
represents more than p% of all the points of the cluster C 2 we project all the arithmetic 
neighbours of x (lying inside or outside the strip t + S 2 ,fc). We choose p such that to 
avoid the superposition of the fully occupied clusters. The projection of a point x with 
n(x) less than p% of all the points of the cluster C 2 is added to the pattern only if it 
is not too close to the already obtained points. We get in this way a discrete set Q 
containing fully occupied copies of the cluster C 2 . 

In the structure analysis of quasicrystals, the experimental diffraction image is 
compared with the diffraction image of the mathematical model, regarded as a set of 
scatterers. In order to compute the diffraction image of a discrete set Q one has to use 
the Fourier transform and to identify il with the Dirac comb ^2 w£Q which will also be 
denoted by Q . The set Q defined in terms of the strip projection method is an infinite 
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For diffraction pattern please see [5] 

C 2 



Figure 3. Left: A set fio containing 923 points defined by starting from a Ci2-clustcr 
C 2 - Centre: The cluster C 2 . Right: For the diffraction pattern of f2 please see [5]. 



For diffraction pattern please see [5] 

C 2 



•' : ':-:{:;(;>v-:}J>::;v:-' : '-fi 

Figure 4. Left: The set fto containing 1019 points corresponding to the set flo from 
the previous figure, defined by using the modified strip projection method. Centre: 
The cluster C 2 - Right: For the diffraction pattern Qq of Qq please see [5]. 



set, but the set Cl we can effectively generate is evidently a finite set, obtained by 
starting from a finite fragment Qq of ^- This is not very bad since any quasicrystal has 
a finite number of atoms. Nevertheless, a fragment of Q or Q can not be an acceptable 
model for a quasicrystal unless it is large enough. 

The diffraction pattern corresponding to f2 is related to the function 



(30) 



where J^[fio] means the Fourier transform of the distribution Q = X^efi I n figure 
3 we present a fragment Qo of the set Q corresponding to a Ci2-cluster and the set 
G K 2 ^o(0 > louo^o(O)} m or d er to illustrate the shape and symmetry properties 
of the diffraction image of f2 . The case of the set f2 corresponding to f2 , obtained by 
using the modified strip projection method in IR 6 , is presented in figure 4. 

It is an open problem if the diffraction properties of the sets obtained by using 
the modified strip projection method are similar to those of sets obtained by the non- 
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modified version. A mathematical answer seems to be difficult, but some suggestions 
in this direction can be obtained by analysing larger fragments. By using our computer 
program the fragment containing 923 points presented in figure 3 and its diffraction 
pattern are obtained in one minute. In two hours one can obtain about 16000 points. 
The fragment generated by the modified method presented in figure 4 can be obtained 
in two minutes. 



4. Quasiperiodic packings of icosahedral clusters 

The icosahedral group Y = 235 can be defined in terms of generators and relations as 
Y = (a,b | a 5 = b 2 = {abf = e) (31) 
and the rotations a, b : 1R 3 — > R 3 

a(a, (3, 7 ) = (*±a - |/3 + ± 7 , |« + \P + ^7, -\<* + ^ + f 7 ) (32) 
6(a,/3,7) = (-a,-/3,7). 

where r = (1 + y / 5)/2, generate an irreducible representation of Y in IR 3 . In the case 
of this representation there are the trivial orbit Y(0, 0, 0) = {(0,0,0)} of length 1, the 
orbits 

Y(a,ar,0) = {g(a,ar, 0) | g G Y} where a G (0, 00) (33) 
of length 12 (vertices of a regular icosahedron) , the orbits 

Y(a,a,a) = {g(a, a, a) \ g G Y} where a G (0, 00) (34) 
of length 20 (vertices of a regular dodecahedron), the orbits 

Y(a, 0, 0) = {g(a, 0, 0) | g G Y} where a G (0, 00) (35) 

of length 30 (vertices of an icosidodecahedron), and all the other orbits are of length 60. 
If the set symmetric with respect to the origin 

C 3 = {vi, v 2 , v k , -v x , -v 2 , -v k } (36) 

where v\ = (i>n, U21, f3i), Vk = (vik, V2k, v^), is a finite union of orbits of Y then the 
vectors 

Wi = (vu,Vi2, Vik) 

w 2 = (v 2 i,v 22 , ...,v 2k ) (37) 
w 3 = (v 3 i,v 32 , ...,v 3k ) 

from M fc are orthogonal 

(w!,w 2 ) = (w 2 ,w 3 ) = (w 3 ,W!) = (38) 

and have the same norm k = \\wi\\ = \\w 2 \\ = \\w 3 \\. They allow us to identify the 
physical space with the three-dimensional subspace 

E 3 = { awi + f3w 2 + 7 w 3 I a, /3, 7 G R } (39) 
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of the superspace The orthogonal projection on E 3 of a vector x G M fc is the vector 
I 3 I = ( I ,^* + ( I *\* + ( I *\* (40) 

The orthogonal projector corresponding to the orthogonal complement 



{ x G R k | (x, y) = for all y G E 3 } 



(41) 



IS 7T3 



R. k — ► E^, tt^x = x — 1T3X. If we describe E3 by using the orthogonal basis 
{k~ 2 wi, k~ 2 W2, k~ 2 w 3 } then the expression in coordinates of 7r 3 is 

7i 3 :R k ^R 3 tt 3 x= {(x,w 1 ),(x,w 2 ),(x,w 3 )). (42) 

The projection W 3j fc = 7r 3 (A^) of the unit hypercube is a polyhedron in the 
(k — 3) -dimensional subspace Eg, and each (k — 4)-dimensional face of W 3fc is the 
projection of a (k — 4)-dimensional face of A*.. Each (k — 4)-face of A^ is parallel to 
(A; — 4) of the vectors ex, e 2 , and orthogonal to four of them. There exist sixteen 

(k — 4)-faces of A k orthogonal to the distinct vectors e^, e i2 , e iz , e i4 , and the set 



X (<^1? *^2> ~^k) 



Xi G {-1/2, 1/2} if i G i 2 , £3, u} 
Xi = if i G" {ii, i 2 , z 3 , i 4 } 

contains one and only one point from each of them. There are 




fc(A; — l)(fc — 2)(fc — 3) 
24 



(43) 



(44) 



sets of 2 4 parallel (k— 4)-faces of A&, and we label them by using the elements of the set 



X 



3.k 



(ii,i2,h,u) G Z 



1 < h < k-3, 
% 2 + 1 < z 3 < k - 1, 



ii + 1 < k < k - 2, 
i 3 + 1 < Z4 < 



(45) 



A point x = (xi, x 2 , Xfc) G M fc belongs to the strip S 3j fc if and only if 



^112*3*4 — 
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H1213M 



(46) 



for each i 2 , i 3 , i±) G X 3) fc, where 



«1*2«3*4 



max 

a,e{-l/2, 1/2} 
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(47) 



For each t G 

p El El EH 



the pattern defined in terms of the strip projection method 



Q = 7r 3 ((t + S 3 , fc ) n Z fc ) = {n 3 x I x-t G S 3ifc and x G Z fc } 
can be regarded as a quasiperiodic packing of copies of the starting cluster C 3 . 



(4f 
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The algorithm based on the strip projection method presented above is very 
efficient. In the case of a three-shell F-cluster formed by the vertices of a regular 
icosahedron, a regular dodecahedron and an icosidodecahedron we use a 31-dimensional 
superspace, W 3 fc is a polyhedron lying in the 28-dimensional subspace Eg bounded 
by 31465 pairs of parallel 2 7- dimensional faces, but we obtain 400-500 points in less 
than 10 minutes [2|. With the modification indicated in the previous section we can 
obtain quasiperiodic packings of multi-shell icosahedral clusters containing a significant 
percentage of fully occupied clusters. 

5. Concluding remarks 

Some of the most remarkable tilings and discrete quasiperiodic sets used in quasicrystal 
physics are obtained by using strip projection method in a superspace of dimension four, 
five or six, and the projection of a unit hypercube as a window of selection |H El Ej- 
The mathematical results presented above allow one to use this very elegant method 
in superspaces of dimension much higher, and to generate discrete quasiperiodic sets 
with a more complicated structure by starting from the symmetry group G and the 
local structure described by a covering cluster C. In our approach the window (which, 
generally, is a polyhedron with hundreds or thousands faces) is described in a simple 
way and we have to compute only determinants of order three or four, independently of 
the dimension of the superspace we use. These mathematical results have allowed us to 
obtain some very efficient computer programs for our algorithm [3] . Hundreds of points 
of our mathematical models can be obtained in only a few minutes. 

The quasiperiodic set generated by starting from a G-cluster C is a packing of 
partially occupied copies of C, but for most of these copies the occupation is very low. 
The main purpose of the paper is to present a modified version of the strip projection 
method. We project certain points lying outside the strip and do not project certain 
points lying inside the strip in order to favour the apperance of fully occupied clusters. 
More exactly, we start from the pattern generated by the standard strip projection 
method and help the clusters with occupation above a certain threshold (p%) to complete 
their configuration up to a fully occupied clusters. We project a minimum number of 
points lying outside the strip, and avoid to project certain points lying inside the strip. 

The symmetry group G corresponding to a real quasicrystal can be deduced from 
the diffraction images, and the covering cluster C can be chosen by analyzing the real 
structure information obtained by high-resolution transmission electron microscopy. 
Our modified strip projection method allows one to generate a mathematical model, 
and to compute the corresponding diffraction image by using the Fourier transform. If 
the agreement with the experimental data is not acceptable one has to look for a more 
suitable covering cluster C. 



Modified strip projection method 



The computer program in FORTRAN 90 used in the case of figure 3 

! PLEASE INDICATE HOW MANY POINTS DO YOU WANT TO ANALYSE 
INTEGER, PARAMETER : : N = 6000 

! PLEASE INDICATE THE DIMENSION M OF THE SUPERSPACE 
INTEGER, PARAMETER : : M = 6 

! PLEASE INDICATE THE RADIUS OF THE PATTERN 
REAL , PARAMETER : : R = 9 . 

INTEGER I, J, K, L, II, 12, 13, JJ, Jl, J2 
REAL Dl, D2, D3, PR 

REAL, DIMENSION (M) :: V, W, TRANSLATION, WJ, EPSILON 
INTEGER, DIMENSION (N) : : CLUSTER 
REAL, DIMENSION (2, M) :: BASIS 
REAL, DIMENSIONS, 2) :: C12 
REAL, DIMENSION^. :M-2, 2 :M-1, 3 :M) :: STRIP 
REAL, DIMENSION (N,M) :: POINTS, STRIPOINTS 
REAL, DIMENSION (N + M) : : XPOINT, YPOINT 
REAL, DIMENSION (200, 200) :: FOURIER 
REAL, DIMENSION (2, 40000) :: PATTERN 
COMPLEX II 
II=(0,1) 

EPSILON = 0.0001 

! PLEASE INDICATE THE COORDINATES OF A POINT OF THE CLUSTER 
BASIS(1,1) = 1.0 
BASIS(2,1) = 0.0 

! PLEASE INDICATE THE TRANSLATION OF THE STRIP YOU WANT TO USE 
TRANSLATION =0.1 

C12(l,l) = SQRT(3.0) / 2.0 
C12(l,2) = -1.0 / 2.0 
C12(2,l) = 1.0 / 2.0 
C12(2,2) = SQRT(3.0) / 2.0 
DO J = 2, 6 
DO I = 1, 2 

BASIS(I,J) = C12(I,1) * BASIS(1,J-1) & 
+ C12(I,2) * BASIS(2,J-1) 
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END DO 
END DO 
STRIP=0 
DO II =1, M-2 
DO 12 =11+1, M-l 
DO 13 =12+1, M 
DO Dl =-0.5, 0.5 
DO D2 =-0.5, 0.5 
DO D3 =-0.5, 0.5 

PR = Dl * BASIS(1,I2) * BASIS(2,I3) + & 
D3 * BASIS(l.Il) * BASIS(2,I2) + & 
D2 * BASISC1.I3) * BASIS(2,I1) - & 
D3 * BASIS(1,I2) * BASISC2.il) - & 

Dl * BASISC1.I3) * BASIS(2,I2) - & 
D2 * BASIS(1,I1) * BASIS(2,I3) 

IF ( PR > STRIP(I1,I2,I3) ) STRIP(I1,I2,I3) = PR 
END DO 
END DO 
END DO 

IF( STRIP(I1,I2,I3) . EQ. ) STRIP (II , 12, I3)=N * SUM( BASIS(1,:) ** 2) 
END DO 
END DO 
END DO 

PRINT*, 'COORDINATES OF THE POINTS OF THE ONE-SHELL C12-CLUSTER: ' 
DO J = 1, M 

PRINT*, J, BASIS(l.J), BASIS(2,J) 
END DO 

PRINT* , ' * STRIP TRANSLATED BY THE VECTOR WITH COORDINATES : ' 
PRINT* , TRANSLATION 

PRINT*, '* PLEASE WAIT A FEW MINUTES OR MORE,& 

DEPENDING ON THE NUMBER OF ANALYSED POINTS' 
POINTS = 
STRIPOINTS = 

POINTSCl,:) = ANINT( TRANSLATION) 
STRIPOINTS ( 1, :) = ANINT( TRANSLATION) 
K = 1 
L = 

DO I = 1, N 

V = POINTSCl, : ) - TRANSLATION 
JJ = 1 

DO II =1, M-2 
DO 12 =11+1, M-l 
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DO 13 =12+1, M 

PR = V(I1) * BASIS(1,I2) * BASIS(2,I3) + & 
V(I3) * BASIS(1,I1) * BASIS(2,I2) + & 

V(I2) * BASIS(1,I3) * BASIS(2,I1) - & 
V(I3) * BASIS(1,I2) * BASIS(2,I1) - & 

V(I1) * BASIS(1,I3) * BASIS(2,I2) - & 
V(I2) * BASIS(1,I1) * BASIS(2,I3) 

IF ( ABS(PR) > STRIP(I1,I2,I3) ) JJ = 

END DO 
END DO 
END DO 

IF( JJ .EQ . 1 ) THEN 
13 = 1 
DO J = 1, L 

WJ = ABS (POINTS (I, :) - STRIPOINTS ( J , : ) ) 
IF( ALL(WJ < EPSILON) ) 13 = 
END DO 

IF( 13 == 1 .AND. SUM( V * V) < R **2 ) THEN 
L = L + 1 

STRIPOINTSCL, :) =POINTS(I,:) 

ELSE 
END IF 

DO II = 1, M 
DO 12 = -1, 1, 2 
W = POINTS (I, : ) 
W(I1) = W(I1) + 12 
13 = 
DO J = 1, K 
WJ = ABS(W - POINTS ( J, :)) 
IF( ALU WJ < EPSILON)) 13 = 1 
END DO 

IF ( 13 == .AND. K < N ) THEN 
K = K + 1 
POINTS (K, : ) = W 
ELSE 
END IF 
END DO 
END DO 
ELSE 
END IF 
END DO 
CLUSTER = 
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DO I = 1, L 
DO Jl = 1, M 

DO J2 = -1, 1, 2 
W = STRIPOINTS(I, :) - TRANSLATION 

W(J1) = W(J1) + J2 
JJ = 1 

DO II =1, M-2 
DO 12 =11+1, M-l 
DO 13 =12+1, M 

PR = W(I1) * BASIS(1,I2) * BASIS(2,I3) + & 
W(I3) * BASIS(l.Il) * BASIS(2,I2) + & 

W(I2) * BASIS(1,I3) * BASIS(2,I1) - & 
W(I3) * BASIS(1,I2) * BASIS(2,I1) - & 

W(I1) * BASIS(1,I3) * BASIS(2,I2) - & 
W(I2) * BASIS(1,I1) * BASIS(2,I3) 

IF ( ABS(PR) > STRIP(I1,I2,I3) ) JJ = 

END DO 
END DO 
END DO 

IF( JJ .EQ . 1 ) CLUSTER(I) = CLUSTER(I) + 1 
END DO 

END DO 

END DO 
DO J = 1, L 

XPOINT(J) = SUM( STRIPOINTSCJ, :) * BASIS(1,:) ) 

YPOINT(J) = SUM( STRIPOINTSCJ,:) * BASIS (2,:) ) 
END DO 

PRINT*, ' NUMBER OF ANALYSED POINTS :', K 

PRINT*, ' NUMBER OF OBTAINED POINTS : ' , L 
DO I = 1, 100 
DO J = 1, 100 
D1=0.0 

DO II = 1, L 

D1=D1+EXP( II * (-1. 5+1*0. 03)*XP0INT(I1)+ & 
II * (-1.5+J*0.03)*YP0INT(I1) ) 

END DO 

F0URIER(I,J)=(ABS(D1))**2 
END DO 
END DO 
12=0 

DO I = 1, 100 
DO J = 1, 100 
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IF(FOURIER(I, J)>0.001*L**2) THEN 
12=12+1 

PATTERN ( 1 , 12) =-1 . 5+1*0 . 03 
PATTERN (2 , I2)=-l . 5+J*0 . 03 
ELSE 
END IF 

END DO 

END DO 

PRINT*, 'INDICATE THE NAME OF A FILE WITH EXTENSION tex FOR RESULTS' 
WRITE (4, 60) 

60 FORMAT ( ' \documentclass{article} & 
\begin{document} & 
\begin{f igure} & 

\setlength{\unitlength}{l . 5mm} & 

\begin{picture}(50,20) (0,0) ' & 
' \put (32 . , 20 . 0) {\circle*{0 .2}} ' ) 
DO J = 1, L 

IF( CLUSTER(J) < M+l ) THEN 
WRITE(4,65) 10+XP0INT(J) , 20+YP0INT(J) 
65 FORMAT ( '\put( ' F10 . 5 ' , ' F10 . 5 , ' ) {\circle*{0 . 2}} ') 
ELSE 
END IF 
END DO 

DO J = 1, L 
IF( CLUSTER(J) > M ) THEN 
WRITE(4,70) 10+XPOINTQ) , 20+YP0INT(J) 
70 FORMAT ( '\put( ' F10 . 5 ' , ' F10 . 5 , ' ) {\circle{0 . 4}} ') 
ELSE 
END IF 
END DO 

DO J = 1, 6 
WRITE(4,72) 32+BASIS(l, J) , 20+BASIS (2 , J) 

72 FORMAT ( '\put( ' F10 . 5 ' , ' F10 . 5 , ' ) {\circle*{0 . 2}} ') 
WRITE(4,73) 32-BASISQ, J) , 20-BASIS (2 , J) 

73 FORMAT ( '\put( ' F10 . 5 ' , ' F10 . 5 , ' ) {\circle*{0 . 2}} ') 
END DO 

WRITE (4, 75) 

75 FORMAT (' \setlength{\unitlength}{l . 8mm} ' ) 
DO J = 1, 12 

WRITE (4, 80) 45+10*PATTERN ( 1 , J) , 17+10*PATTERN (2 , J) 
80 FORMAT ( '\put( ' F10 . 5 ' , ' F10 . 5 , ' ) {\circle*{0 . 1}} ') 
END DO 
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WRITE (4, 90) 
90 FORMAT ( ' \end{picture} & 

\caption-[Quasiperiodic set obtained by using & 

the strip projection method } & 
\end{f igure} & 
\end{document} ' ) 
PRINT*, '* COMPILE THE OBTAINED FILE AND SEE THE ".dvi" FILE' 
END 
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The computer program in FORTRAN 90 used in the case of figure 4 

! PLEASE INDICATE HOW MANY POINTS DO YOU WANT TO ANALYSE 
INTEGER, PARAMETER : : N = 6000 

! PLEASE INDICATE THE DIMENSION M OF THE SUPERSPACE 
INTEGER, PARAMETER : : M = 6 

! PLEASE INDICATE THE RADIUS OF THE PATTERN 
REAL , PARAMETER : : R = 9 . 

INTEGER I, J, K, L, II, 12, 13, JJ, Jl, J2, LI 
REAL Dl, D2, D3, PR, XP, YP 

REAL, DIMENSION (M) :: V, W, TRANSLATION, WJ, EPSILON 
INTEGER, DIMENSION (N) : : CLUSTER 
REAL, DIMENSION (2, M) :: BASIS 
REAL, DIMENSIONS, 2) :: C12 
REAL, DIMENSION^. :M-2, 2 :M-1, 3 :M) :: STRIP 
REAL, DIMENSION (N,M) :: POINTS, STRIPOINTS 
REAL, DIMENSION (N + M) : : XPOINT, YPOINT 
REAL, DIMENSION (200, 200) :: FOURIER 
REAL, DIMENSION (2, 40000) :: PATTERN 
COMPLEX II 
II=(0,1) 

EPSILON = 0.0001 

! PLEASE INDICATE THE COORDINATES OF A POINT OF THE CLUSTER 
BASIS(1,1) = 1.0 
BASIS(2,1) = 0.0 

! PLEASE INDICATE THE TRANSLATION OF THE STRIP YOU WANT TO USE 
TRANSLATION =0.1 

C12(l,l) = SQRT(3.0) / 2.0 
C12(l,2) = -1.0 / 2.0 
C12(2,l) = 1.0 / 2.0 
C12(2,2) = SQRT(3.0) / 2.0 
DO J = 2, 6 
DO I = 1, 2 

BASIS(I,J) = C12(I,1) * BASIS(1,J-1) & 
+ C12(I,2) * BASIS(2,J-1) 
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! BASIS(I,6+J) = C12(I,1) * BASIS(1,5+J) & 
! + C12CI.2) * BASIS(2,5+J) 

END DO 
END DO 
STRIP=0 
DO II =1, M-2 
DO 12 =11+1, M-l 
DO 13 =12+1, M 
DO Dl =-0.5, 0.5 
DO D2 =-0.5, 0.5 
DO D3 =-0.5, 0.5 

PR = Dl * BASIS(1,I2) * BASIS(2,I3) + & 
D3 * BASIS(l.Il) * BASIS(2,I2) + & 
D2 * BASIS(1,I3) * BASIS(2,I1) - & 
D3 * BASIS(1,I2) * BASIS(2,I1) - & 

Dl * BASIS(1,I3) * BASIS(2,I2) - & 
D2 * BASIS(1,I1) * BASIS(2,I3) 

IF ( PR > STRIP(I1,I2,I3) ) STRIP(I1,I2,I3) = PR 
END DO 
END DO 
END DO 

IF( STRIP(I1,I2,I3) . EQ. ) STRIP(I1 , 12, I3)=N * SUM( BASIS(1,:) ** 2) 
END DO 
END DO 
END DO 

PRINT*, ' COORDINATES OF THE POINTS OF THE ONE-SHELL C12-CLUSTER : ' 
DO J = 1, M 

PRINT*, J, BASIS(l.J), BASIS(2,J) 
END DO 

PRINT* , ' * STRIP TRANSLATED BY THE VECTOR WITH COORDINATES : ' 
PRINT*, TRANSLATION 

PRINT*, '* PLEASE WAIT A FEW MINUTES OR MORE,& 

DEPENDING ON THE NUMBER OF ANALYSED POINTS' 
POINTS = 
STRIPOINTS = 

POINTSCl,:) = ANINT( TRANSLATION) 
STRIPOINTS ( 1, :) = ANINT( TRANSLATION) 
K = 1 
L = 

DO I = 1, N 

V = POINTSCl, : ) - TRANSLATION 
JJ = 1 
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DO II =1, M-2 
DO 12 =11+1, M-l 
DO 13 =12+1, M 

PR = V(I1) * BASIS(1,I2) * BASIS(2,I3) + & 
V(I3) * BASIS(l.Il) * BASIS(2,I2) + & 

V(I2) * BASIS(1,I3) * BASIS(2,I1) - & 
V(I3) * BASIS(1,I2) * BASIS(2,I1) - & 

V(I1) * BASIS(1,I3) * BASIS(2,I2) - & 
V(I2) * BASIS(1,I1) * BASIS(2,I3) 

IF ( ABS(PR) > STRIP(I1,I2,I3) ) JJ = 

END DO 
END DO 
END DO 

IF( JJ .EQ. 1 ) THEN 
13 = 1 
DO J = 1, L 

WJ = ABS (POINTS (I, :) - STRIPOINTS ( J , : ) ) 
IF( ALL(WJ < EPSILON) ) 13 = 
END DO 

IF( 13 == 1 .AND. SUM( V * V) < R **2 ) THEN 
L = L + 1 

STRIPOINTS (L, :) = POINTS (I,:) 

ELSE 
END IF 

DO II = 1, M 
DO 12 = -1, 1, 2 
W = POINTS (I, : ) 
W(I1) = W(I1) + 12 
13 = 
DO J = 1, K 
WJ = ABS(W - POINTS ( J, :)) 
IF( ALU WJ < EPSILON)) 13 = 1 
END DO 

IF ( 13 == .AND. K < N ) THEN 

K = K + 1 

POINTS (K, : ) = W 

ELSE 

END IF 
END DO 
END DO 
ELSE 
END IF 
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END DO 
CLUSTER = 
DO I = 1, L 
DO Jl = 1, M 

DO J2 = -1, 1, 2 
W = STRIPOINTSCI, :) - TRANSLATION 

W(J1) = W(J1) + J2 
JJ = 1 

DO II =1, M-2 
DO 12 =11+1, M-l 
DO 13 =12+1, M 

PR = W(I1) * BASIS(1,I2) * BASIS(2,I3) + & 
W(I3) * BASIS(l.Il) * BASIS(2,I2) + & 

W(I2) * BASIS(1,I3) * BASIS(2,I1) - & 
W(I3) * BASIS(1,I2) * BASIS(2,I1) - & 

W(I1) * BASIS(1,I3) * BASIS(2,I2) - & 
W(I2) * BASIS(1,I1) * BASIS(2,I3) 

IF ( ABS(PR) > STRIP(I1,I2,I3) ) JJ = 

END DO 
END DO 
END DO 

IF( JJ .EQ . 1 ) CLUSTER(I) = CLUSTER(I) + 1 
END DO 

END DO 

END DO 
II = L 

DO I = 1, II 

IF ( CLUSTER(I) > M) THEN 
DO Jl = 1, M 
DO J2 = -1, 1, 2 
W = STRIPOINTSCI, :) 
W(J1) = W(J1) + J2 
13 = 
DO J = 1, L 
WJ = ABS(W - STRIPOINTSCJ, :)) 
IF( ALL( WJ < EPSILON)) 13 = 1 
END DO 
IF ( 13 == ) THEN 
L = L + 1 

STRIPOINTSCL, : ) = W 

ELSE 

END IF 
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END DO 
END DO 
ELSE 
END IF 
END DO 
L1=0 
DO I = 1, II 

IF ( CLUSTER(I) > M) THEN 
L1=L1+1 

XPOINT(Ll) = SUM( STRIPOINTSU, :) * BASIS(1,:) ) 

YPOINT(Ll) = SUM( STRIPOINTSCI,:) * BASIS(2,:) ) 

ELSE 
END IF 
END DO 

DO J = 11+1, L 
L1=L1+1 

XPOINT(Ll) = SUM( STRIPOINTSCJ,:) * BASIS(1,:) ) 

YPOINT(Ll) = SUM( STRIPOINTSCJ, : ) * BASIS(2,:) ) 
END DO 
Dl=4.0 
DO I =2, M 

IF(((BASIS(1,1)+BASIS(1,I))**2 + (BASIS(2, 1)+BASIS(2,I))**2 ) < Dl ) & 
D1=(BASIS(1,1)+BASIS(1,I))**2 + (BASIS (2 , 1) +BASIS (2 , I) ) **2 
END DO 
DO I =2, M 

IF(((BASIS(1,1)-BASIS(1,I))**2 + (BASIS (2 , 1) -BASIS (2 , I) ) **2 ) < Dl ) & 
D1=(BASIS(1,1)-BASIS(1,I))**2 + (BASIS (2 , 1) -BASIS (2 , I) ) **2 
END DO 

DO I = 1, II 

IF ( CLUSTER(I) <= M) THEN 
D2=4 . 

XP = SUM( STRIPOINTSCI,:) * BASIS(1,:) ) 

YP = SUM( STRIPOINTSCI,:) * BASIS (2,:) ) 
DO J = 1, LI 

IFC ((XP - XP0INT(J))**2 + (YP - YP0INT(J))**2 ) < D2 ) & 
D2=(XP - XP0INT(J))**2 + (YP - YPOINT(J) ) **2 
END DO 

IFC D2 > 0.9*D1 ) THEN 

L1=L1+1 

XPOINT(Ll)=XP 

YPOINT(Ll)=YP 

ELSE 
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END IF 

ELSE 
END IF 
END DO 

DO I = 1, 100 
DO J = 1, 100 
D1=0.0 

DO II = 1, LI 

D1=D1+EXP( II * (-1. 5+1*0. 03)*XP0INT(I1)+ & 
II * (-1.5+J*0.03)*YP0INT(I1) ) 

END DO 

FOURIER (I, J) =(ABS(D1))**2 
END DO 
END DO 
12=0 

DO I = 1, 100 
DO J = 1, 100 

IF(FOURIER(I , J) >0 . 0015*L**2) THEN 
12=12+1 

PATTERN ( 1 , 12) =-1 . 5+1*0 . 03 
PATTERN (2 , 12) =-1 . 5+ J*0 . 03 
ELSE 
END IF 

END DO 

END DO 

PRINT*, ' NUMBER OF ANALYSED POINTS :', K 

PRINT*, ' NUMBER OF OBTAINED POINTS : ' , LI 

PRINT*, 'INDICATE THE NAME OF A FILE WITH EXTENSION tex FOR RESULTS 
WRITE (4, 60) 

60 FORMAT C\documentclass{article} & 
\begin{document} & 
\begin{f igure} & 

\setlength{\unitlength}{l . 5mm} & 

\begin{picture}(50,20) (0,0) ' & 
' \put (32 . , 20 . 0) {\circle*{0 .2}} ' ) 
DO J = 1, LI 

IF( CLUSTER(J) < M+l ) THEN 
WRITE(4,65) 10+XPOINTQ) , 20+YP0INT(J) 
65 FORMAT ( '\put( ' F10 . 5 ' , ' F10 . 5 , ' ) {\circle*{0 . 2}} ') 
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ELSE 
END IF 
END DO 

DO J = 1, LI 
IF( CLUSTER(J) > M ) THEN 
WRITE(4,70) 10+XPOINTO) , 20+YP0INT(J) 
70 FORMAT ( '\put( ' F10 . 5 ' , ' F10 . 5 , ' ) {\circle*{0 . 2}} ') 
ELSE 
END IF 
END DO 

DO J = 1, 6 
WRITE(4,72) 32+BASISC1, J) , 20+BASIS (2 , J) 

72 FORMAT ( '\put( ' F10 . 5 ' , ' F10 . 5 , ' ) {\circle*{0 . 2}} ') 
WRITE(4,73) 32-BASISU, J) , 20-BASIS (2 , J) 

73 FORMAT ( '\put( ' F10 . 5 ' , ' F10 . 5 , ' ) {\circle*{0 . 2}} ') 
END DO 

WRITE (4, 75) 

75 FORMAT (' \setlength{\unitlength}{l . 8mm} ' ) 
DO J = 1, 12 

WRITE (4, 80) 45+10*PATTERN(l, J) , 17+10*PATTERN (2 , J) 
80 FORMAT ( '\put( ' F10 . 5 ' , ' F10 . 5 , ' ) {\circle*{0 . 1}} ') 
END DO 

WRITE(4,90) 
90 FORMAT ( ' \end{picture} & 

\caption-COuasiperiodic set obtained by using & 
the modified strip projection method } & 
\end{f igure} & 
\end{document} ' ) 
PRINT*, '* COMPILE THE OBTAINED FILE AND SEE THE ".dvi" FILE' 
END 
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